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Abstract. Planet accretion onto solar type stars may significantly change the stellar atmospheric abundances of 
^Li and ^Li if it takes place after the star has arrived at the main sequence. Ingestion of planets at an earlier 
phase will not affect theses abundances because of extensive pre-main sequence mixing and burning. We present 
quantitative estimates of the main sequence evolution of stellar surface lithium abundances after planet ingestion. 
At solar metallicities, for stellar masses in the range ~ 1.3 — 1.1 M0, a large fraction of the ingested ®Li is likely 
to be preserved during the whole main sequence lifetime according to the standard model. Preservation of the 
less fragile ^Li isotope occurs in a larger mass range ~ 1.3 — 0.9 M©. At high metallicities typical of planet host 
stars, the ranges of masses are slightly reduced essentially due to the thicker convective zones. Further reduction 
is expected if non-standard processes cause extra-mixing of material below the base of the convective zone, but 
even in this case there are stellar masses where both isotopes are significantly preserved. We conclude that there 
is a range of stellar effective temperature where migration and accretion of planets and planetary material can be 
empirically tested using high-resolution spectroscopy of the lithium isotopes. 
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""i 1. Introduction 



The accretion of planets or planetary material onto stars 
have been considered as a possible explanation of the high 
metallicity frequently found in planet host stars (Laughlin 
& Adams 1997; Gonzalez 1997, 1998, 2001; Sandquist et 
al. 1998; Laughlin 2000; Santos et al. 2000; Murray et 
al. 2001; Pinsonneault et al. 2001). The very recent de- 
tection of the ^Li isotope in the atmosphere of the solar- 
type star HD 82943 (Israelian et al. 2001), which hosts a 
planetary system (Mayor et al. 2001), has also been in- 
terpreted as evidence for accretion of planets. Jupiter-like 
planets are expected to preserve their original content of 
the lithium isotopes, while stars like HD 82943 with a 
mass ^ 1.1 Mq destroy, via nuclear reactions with pro- 
tons, all its ^Li (and a small fraction of the more robust 
'^Li) during the pre-main sequence (PMS) evolution (see 
e.g. Proffitt and Michaud 1989). It has been argued that 
ingestion and dilution of planetary material in HD 82943 
after its superficial convective zone receded to the main- 
sequence location allowed long-term preservation of the 
newly acquired abundances of lithium isotopes. In fact, 
according to standard models for this type of stars, the 



base of the convective zone does not reach the tempera- 
tures needed for ^Li and ^Li burning (T > 2.2 x 10^ K and 
T > 2.6 X 10^ K, respectively). The rate of destruction of 
^Li by protons has a temperature dependence similar to 
that of ^Li, but on much shorter time scale (i.e. 100 times 
shorter at 1 Mq in the standard model). 

In this letter, we examine in quantitative terms the 
possibility that the dilution of one or several Jupiter-like 
planets in the convective envelope of a main sequence 
star could appreciably change the observed abundance of 
^Li and ^Li . In the framework of the standard model we 
predict the evolution of these isotopic abundances in stars 
with masses ranging from 1.3 to 0.6 Mq. We also discuss 
the effect due to some transport processes which cause 
extra-mixing below the convective zone (see for review 
Michaud & Charbonneau 1991). We will show that there 
is a narrow range of mass (or Tog) where the accretion of 
planetary material could produce an observable increase 
of ^'^Li abundance during a significant fraction of the life- 
time of the star. 
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2. Models 

The stellar models have been computed with the up- 
dated version of the ATON2.0 code (Ventura et al. 1998). 
We constructed models with masses ranging from 0.6 to 
1.3 Mq with a step of AM = 0.05 M©, and followed the 
evolution from the pre-main sequence to 7 Gyr. These 
models use a mixing length theory (MLT) description of 
convective transport with a parameter of mixing length 
a = 1.6, as given by the solar calibration (Ventura et 
al. 1998; Montesinos et al. 2001). Two different chem- 
ical compositions were adopted: a) initial solar helium 
abundance of Yq — 0.28 and heavy-elements mass frac- 
tion of Zq = 0.02 ([Fe/H]=0.0) (Montesinos et al. 2001), 
h) Z ^ 0.04, with AY/AZ = 2, that is, [Fe/H]=-|-0.34. 

The code used in our model computations allows to 
include gravitational settling of helium and thermal dif- 
fusion, but it does not consider the diffusion of metals. 
Diffusion generally produces a reduction of the convec- 
tion zone, however, Turcotte et al. (1998) showed that 
the convective envelopes in models with metal diffusion 
are deeper than in models including only helium settling. 
Furthermore, their results indicate (see their Fig. 4) that 
for the higher masses considered in this paper (Af < 
1.3 Mq) the mass of the convective envelope of a non- 
diffusion model is actually closer to the full-diffusion one 
than to a model incorporating He settling only (in fact, the 
non-diffusion and full-diffusion curves overlap) . Because of 
that, we decided to use stellar structures obtained with- 
out including diffusion. Notice also that, a proper treat- 
ment of diffusion (Turcotte et al. 1998) has no significant 
impact on the evolution of lithium abundance in stars 
with M < 1.3 Mq. The nuclear burning rates of ^Li and 
^Li in the star's convective envelope were computed for 
each temporal step in evolution using NACRE nuclear re- 
action tables (Angulo et al. 1999). 

Masses for the engulfed sub-stellar companions ranged, 
from 1 Mj0 to 40 Mj . This upper mass limit was adopted 
to ensure the complete preservation of ^Li (Pozio 1991, 
Nelson et al. 1993), but, in practice, the maximum mass 
of the accreted planet was constrained to be less than the 
mass in the convective zone of the star under considera- 
tion. 

The hypothesis and approximations adopted in our 
study are: i) planet accretion takes place just after the star 
arrived on the main sequence. The analysis presented here 
does not change if the planet is engulfed at any other mo- 
ment during the main sequence lifetime. Such case can be 
dealt by taking into account that the temperature, mass 
of the convective envelope and the stellar surface lithium 
abundance will be different since they decrease as age in- 
creases, ii) The planet is completely dissolved in the stellar 
convective envelope. As shown by numerical simulations 
(Sandquist et al. 1998), the reliability of this assumption 
depends on the mass of the planet (Mpi) and on the mass 
of the convective envelope Mcz, and also strongly depends 

^ Jupiter mass, Mj = 9.55 x lO^^M© 



on the details of the internal structure of the planet. The 
smaller the Mqz , the lower the fraction of planet dissolved 
in the convective zone. That means that, at least for the 
higher planet masses under consideration, the quantity of 
lithium provided by the planet must be reduced in a factor 
equal to the fraction of the planet consumed in the star's 
convective envelope. Hi) The engulfment of a planet has 
no significant effect in the structure of the star, iv) The 
planet is instantly mixed in the convective zone of the star. 
Given that the diffusion coefficient in the convective zone 
is of the order of lO^^cm^s"^, the mixing time is lower 
than the temporal step in the evolution computation. 

The stellar models were followed from the PMS, but 
we will not use the surface abundance of ^Li and ^Li given 
by the models because it is well known that they overes- 
timate PMS lithium destruction (D'Antona & Mazzitelli 
1997). For masses lower than 1.4 Mq it is expected that 
no ^Li exists after the PMS (Proflltt & Michaud 1989). 
As regards surface ^Li abundance we adopted the values 
observed in stars of the Pleiades Cluster (Soderblom et 
al. 1993). So the stellar surface abundances of lithium iso- 
topes at the moment of the engulfment of the planet are: 

iV(6Li)/A^(iJ) = and NCU)/N{H) = i^(reff)pieiades. 

These abundances in the planets are assumed equal to 
the meteoritic values. That is: NCU)o/N{H) = 2.10"^, 
and iV(^Li)o/iV(^Li)o = 12. 

Following Burrows & Sharp (1999) we assume an 
Anders & Grevesse (1989) solar chemical composition for 
both the planet and the star. In this case, the difference 
between the abundance of the isotope J in the star's con- 
vective zone before (s) and after (s+p) the dilution of the 
planet can be written as: 

e( J).+p - e( J). = log (^-^ •(!-/))- log / . (1) 

Where e(J) ^ log{N (J) / N {H)) -t- 12, a; = Afpi/Mcz, and 
/ is the depletion factor of J in the convective zone just 
before the absorption of the planet. So the change in the 
stellar surface abundance depends strongly on the factor 
/, and it is evident from Eq. (|l|) that the lithium abun- 
dance enhancement will be significant only if the mass 
fraction of the isotope J in the planet is much higher than 
in the stellar convective envelope. If the engulfed planet 
has lost a significant part of its initial content of Hydrogen, 
the mass fraction of lithium in the planet is larger than in 
the solar chemical composition case. Therefore, also the 
enhancement of atmospheric lithium abundance will be 
larger, and the results presented here would be a lower 
limit. 

3. Evolution of Li provided by planetary material 

3.1. Standard model framework 

The standard model assumes that no exchange of mate- 
rial occurs between the convective envelope and the stable 
radiative interior of the star. Therefore, the surface chem- 
ical composition reflects only the composition of the well 
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Fig. 1. ^Li (lower curves) and ^Li (upper curves) life- 
time against (p, a) nuclear reactions [N [Li) j N (Li)q = 
exp(— i/rw)) as a function of effective temperature. The 
curves correspond to standard models at 10* yr with 
two different metallicities: [Fe/H]=0.0 (solid lines) and 
[Fe/H]=-|-0.34 (dashed lines). The numbers close to the 
symbols (empty circles for solar metallicity and full circles 
for [Fe/H]=-|-0.34) indicate the mass of the stellar model. 



mixed convective envelope. In this subsection we present 
the evolution of lithium in the stellar surface due to the 
nuclear burning taking place in the convective envelope. 

Assuming solar chemical composition, a planet of 1 Mj 
supplies ^ 1.4 10'*'' atoms of ^Li and ~ 1.7 10^^ atoms of 
''Li. What is the evolution of these atoms during the main 
sequence lifetime of the star? The nuclear time scale, tn , is 
an estimation of how long these atoms could survive to the 
nuclear burning in the convective envelope, tn, is defined 
as the time needed to decrease the superficial abundance 
of the the isotope J by a factor e, and it has been ob- 
tained from the nuclear reaction rates at each layer in the 
convective zone, and assuming instantaneous mixing. In 
Fig. |l| we plot Tn for ^Li and ''Li at 100 Myr as a function 
of effective temperature. We plot the curves correspond- 
ing to two different values of metallicity. We see that for 
logToff > 3.76 (~ 5700 K), a large fraction of the ^Li from 
the planet will survive for a long time in the convective 
envelope of the star. The value of tn given in the plot 
is a minimum, since during the main sequence lifetime, 
the temperature at the bottom of the convective zone de- 
creases, and so the nuclear time scale increases. The effect 
of increasing the metallicity is to reduce the nuclear time 
scale at a given Teff. A factor 2 higher metallicity implies 
a factor 3 lower time scale at 5500 K, a factor 2 at 5700 K 
and only a reduction of 10% for Toff higher than 5800 K. 



Fig. II presents the evolution (based on models and 
hypothesis described in Sect. 2) of ^Li and ^Li abundance 
for several stellar masses, and for two values of metal- 
licity: [Fe/H]=0.0 (left panel) and [Fe/H]=+0.34 (right 
panel) . That temporal dependence is obtained computing 
the nuclear depletion of ^Li and ^Li in the stellar convec- 
tive envelope for each temporal step in the evolution and 
assuming instantaneous mixing. We plot only values of 
lithium abundances corresponding to times longer than 
the time of planet accretion, and lower than the drag-up 
of the convective envelope. In the upper panels, we plot 
the curves corresponding to the ^Li that result from the 
accretion of a planet with ~ 1, 5, and 10 Mj in stars with 
masses between 0.9 and 1.2 Mq (for 1.25 and 1.3 Mq the 
enhancement of ^Li is slightly larger because the mass of 
the convective envelope is smaller, and the curves are par- 
allel to those of 1.2 M0). This preliminary and simple 
modelling provides strong support for the interpretation 
by IsraeUan et al. (2001) of the ^Li detection in HD 82943. 
According to our models, assuming a Tcff ^ 6000 K, age 
of 5 Gyr and metaUicity [Fe/H]- 0.32 (Santos et al. 2000) 
this star has a mass of ~ 1.2 M©. As can be seen in Figg 
right panel, it should preserve any ingested ^Li during all 
the main sequence lifetime. 

The ^Li abundances plotted in Fig. ^ were obtained 
assuming a complete dissolution of the planet in the 
stellar convective zone. For solar metallicity, the masses 
of the convective envelope in stars with 1.2, 1.1 and 

1.0 Mq are 0.006 Mq, 0.017 Mq and 0.034 Mq, respec- 
tively. Numerical simulations by Sandquist et al. (1998) 
indicate that as the mass of the convective envelope de- 
creases, the fraction of planet that is dissolved is reduced. 
If only 0.3 of a Jupiter-like planet is dissolved in the con- 
vective zone of a 1.2 star, the corresponding curve of 
^Li will only decrease in ^ 0.2 dex. 

Fig. H shows that there is a mass range in which 
^Li could survive. In fact, for stars with masses larger than 

1.1 Mq (for Zq, or 1.15 Mq for Z = 2Zq) the superficial 
abundance of ^Li does not change with age. The reason is 
that during their main sequence lifetime, the temperature 
at the base of the convective zone is not high enough to 
produce efficient burning of ®Li. As the stellar mass de- 
creases, the depth of the convective zone increases, and 
then, even if a large quantity of planetary material is dis- 
solved in the stellar envelope, ^Li will be depleted in a 
short time scale ('^ 10* yr). 

The lower panels of Fig. || presents the predictions for 
^Li according to standard model. Note that, opposite to 
^Li, the ^Li isotope is depleted on the PMS only by a small 
factor, and that stars with masses larger than O.OM© will 
not deplete a significant amount during the main sequence 
lifetime. If these stars were able to dissolve 10 Mj in their 
convective envelope, the effect on the abundance would be 
to introduce only a small enhancement, 0.2 dex at most. 
The reason is shown in Eq. (|l|), where we see that the 
increase of abundance depends strongly on the depletion 
factor of ''Li in the convective envelope before the accre- 
tion of the planetary material. 
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Fig. 2. Evolution of ^Li (upper panels) and ''Li (lower panels) provided by dissolution of planets in the convective 
envelope of stars with masses from 1.2 to 0.7 Mq (numbers indicate the mass of the star) and metallicities : [Fe/H]=0.0 
(left panels) and [Fe/H]=+0.34 (right panels). ^Li : for each stellar model there are three curves corresponding to 
(from down to top) ~ 1, 5 and 10 Mj of planetary material dissolved. ^Li : as for ^Li, but we also plot the curve 
corresponding to no planet accretion. The ^Li curves corresponding to : 1.0, 1.1 and 1.2 Mq (left panel), and 1.1, 1.15 
and 1.2 Mq (right panel) are concentrated in the top of the respective figures. 



3.2. Impact of non-standard transport processes 

Observations of lithium in open clusters show that the at- 
mospheric abundance decreases during the main sequence 
lifetime. In order to explain this main-sequence depletion 
different mixing processes are invoked (e.g. Michaud & 
Charbonneau 1991), and we should consider their effect 
on the evolution of ®Li in the stellar mass range where 
standard models predict complete preservation. Among 
these processes are: a) microscopic diffusion, that consists 
roughly in the sink of heavy elements with respect to the 
light ones due to the gravitational potential; b) turbulent 
diffusion which causes mixing of material between the con- 
vective zone and the deep layers where the temperature is 
high enough to burn lithium. 

Concerning the first one, the diffusion velocities are 

similar for ^Li and ^Li, but effects of microscopic diffu- 
sion are negligible for the effective temperature and ages 
considered here (e.g. Richer & Michaud 1993). Turcotte 
et al. (1998) showed that for 1.3 Mq star, the effect of mi- 
croscopic diffusion on Li is at most a depletion by a fac- 



tor 0.55, and only after 2 Gyr. In order to test the effect 
of turbulent mixing we have considered a diffusion coef- 
ficient produced by internal waves that fits the evolution 
of ^Li abundances between the Pleiades and M67 clus- 
ters quite well (Montalban & Schatzman 2000): ^Li would 
decrease only in 0.2 dex after 5 Gyr in a 1.2 Mq star 
(Zq). However, in a 1.1 Mq {Zq) star, ^Li will decrease 
0.6 dex in the first Gyr and 1.7 dex after 5 Gyr. Of 
course, these results depend strongly on the behaviour of 
the diffusion coefficient close to the base of the convec- 
tive zone. We recall that the diffusion coefficient produced 
by internal waves has large amplitude at the CZ bound- 
ary, and then it decreases rapidly downwards. Therefore, 
the results are very sensitive to the distance between the 
base of the convective zone and the ^Li burning layer. We 
have also tested a parametric turbulent diffusion coeffi- 
cient Dt = Do{pcz/p)^ as proposed by Proffitt & Michaud 
(1991) {pcz is the density at the base of convective zone). 
A value of Dq ~ 1.510^cm^s~^ produces, after 5 Gyr, a 
^Li depletion of : 1.6 dex for 1 Mq (~ 5800 K), 0.2 dex 
for 1.1 M0 (~ 6000 K), and < 0.1 dex for 1.2 Mq (~ 
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6100 K) (these values are in good agreement with estima- 
tions by Pinsonneault et al. 2001 using mixing induced 
by rotation), and causes a ^Li depletion of 0.9 dex at 
Teff ~6000 K, 0.05 dex at T^s ~6100 K, and a complete 
depletion after 1 Gyr for T^g ~5800 K. 

From observations it is inferred that the efficiency of 
non-standard mixing decreases with the age in the main- 
sequence. So what happens if the accretion of the plane- 
tary material takes place in a ^Li depleted star at a suf- 
ficently old age for extra-mixing to be inefficient? The 
stellar ''Li abundance could be significanly enhanced for 
a long period of time. For instance, if a 1 Mj planet were 
dissolved in our Sun's convective envelope would produce 
an enhancement of 0.3 dex, if 10 Mj were accreted the 
corresponding increase would be 0.9 dex. 

4. Conclusions 

The results presented here show that the accretion of 
planetary material by a solar-type star will produce an 
enhancement of the stellar surface abundance of lithium 
which could be observable during a large fraction of the 
main sequence lifetime. Within the framework of the stan- 
dard model, the Toff range where the ingested ^Li could be 
detected is 5900-6400 K. The upper limit is imposed by 
observational evidence given by the Li-Gap in the Hyades 
(Boesgaard & Tripicco 1986). The lower limit is the result 
of direct convective mixing and should increase about 100 
or 200 K when extra-mixing processes are considered. The 
Teff range where ^Li would survive is larger. Therefore, a 
very valuable test for planetary migration/accretion sce- 
narios can be provided by high-resolution spectroscopic 
observations of ^•''Li in planet host stars. 
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